Molecular characterization of cetacean poxviruses along the coast of mainland Portugal.

Cetacean poxvirus (CePV) is the causative agent of tattoo skin disease (TSD) in dolphins, porpoises and whales, a condition characterized by pinhole, ring-like lesions or generalized tattoo-like skin lesions. This study genetically characterized cetacean poxviruses from stranded animals along mainland Portugal. Samples from skin lesions compatible with TSD were obtained from 4 odontocete species (Delphinus delphis, Stenella coeruleoalba, Phocoena phocoena, and Tursiops truncatus) and analyzed using a conventional PCR assay targeting the DNA polymerase gene partially. Among the positive samples (n = 29, 65.9%), a larger DNA polymerase gene fragment was obtained, allowing a robust phylogenetic analysis. Nineteen samples (43.2%) were successfully amplified and sequenced using Sanger sequencing. By combining 11 of these sequences with those from public databases, a maximum likelihood phylogenetic tree was constructed, revealing high heterogeneity within the group. These findings contribute to a better understanding of the genetic diversity, epidemiology, phylogenetics, and evolution of CePV.

Human orf virus (family Poxviridae) infection following a lamb bite in Hungary.

Human orf disease (called ecthyma contagiosum or contagious/infectious pustular dermatitis in animals) was confirmed on the fingers of both hands of a 24-year-old female, after feeding diseased lambs with a nursing bottle in April 2023. In addition to skin symptoms, she had low-grade fever (37.6°C) and swollen lymph nodes in both axilla. The presence of orf virus (genus Parapoxvirus, family Poxviridae) was confirmed, and this strain, Baja/2023/HUN (OR372161-OR372163), was found to have > 98% nucleotide sequence identity to sheep-origin orf viruses in four tested genome regions (ORF011/B2L, ORF019, ORF020/VIR, and ORF056). This is the first report of a human case of infection with the neglected zoonotic orf virus in Hungary.

Poxvirus A51R proteins regulate microtubule stability and antagonize a cell-intrinsic antiviral response.

Numerous viruses alter host microtubule (MT) networks during infection, but how and why they induce these changes is unclear in many cases. We show that the vaccinia virus (VV)-encoded A51R protein is a MT-associated protein (MAP) that directly binds MTs and stabilizes them by both promoting their growth and preventing their depolymerization. Furthermore, we demonstrate that A51R-MT interactions are conserved across A51R proteins from multiple poxvirus genera, and highly conserved, positively charged residues in A51R proteins mediate these interactions. Strikingly, we find that viruses encoding MT interaction-deficient A51R proteins fail to suppress a reactive oxygen species (ROS)-dependent antiviral response in macrophages that leads to a block in virion morphogenesis. Moreover, A51R-MT interactions are required for VV virulence in mice. Collectively, our data show that poxviral MAP-MT interactions overcome a cell-intrinsic antiviral ROS response in macrophages that would otherwise block virus morphogenesis and replication in animals.

[Poxvirus-encoded DNA replication proteins: potential targets for antivirals]. / Le complexe de réplication des poxvirus : cible potentielle de molécules antivirales.

In the spring of 2022, an epidemic due to human monkeypox virus (MPXV) of unprecedented magnitude spread across all continents. Although this event was surprising in its suddenness, the resurgence of a virus from the Poxviridae family is not surprising in a world population that has been largely naïve to these viruses since the eradication of the smallpox virus in 1980 and the concomitant cessation of vaccination. Since then, a vaccine and two antiviral compounds have been developed to combat a possible return of smallpox. However, the use of these treatments during the 2022 MPXV epidemic showed certain limitations, indicating the importance of continuing to develop the therapeutic arsenal against these viruses. For several decades, efforts to understand the molecular mechanisms involved in the synthesis of the DNA genome of these viruses have been ongoing. Although many questions remain unanswered up to now, the three-dimensional structures of essential proteins, and in particular of the DNA polymerase holoenzyme in complex with DNA, make it possible to consider the development of a model for poxvirus DNA replication. In addition, these structures are valuable tools for the development of new antivirals targeting viral genome synthesis. This review will first present the molecules approved for the treatment of poxvirus infections, followed by a review of our knowledge of the replication machinery of these viruses. Finally, we will describe how these proteins could be the target of new antiviral compounds.

Tracing the journey of poxviruses: insights from history.

The historical significance of the poxviruses is profound, largely due to the enduring impact left by smallpox virus across many centuries. The elimination of smallpox is a remarkable accomplishment in the history of science and medicine, with centuries of devoted efforts resulting in the development and widespread administration of smallpox vaccines. This review provides insight into the pivotal historical events involving medically significant poxviruses. Understanding the remarkable saga of combatting smallpox is crucial, serving as a guidepost for potential future encounters with poxvirus infections. There is a continual need for vigilant observation of poxvirus evolution and spillover from animals to humans, considering the expansive range of susceptible hosts. The recent occurrence of monkeypox cases in non-endemic countries stands as a stark reminder of the ease with which infections can be disseminated through international travel and trade. This backdrop encourages introspection about our journey and the current status of poxvirus research.

Neutralization Determinants on Poxviruses.

Smallpox was a highly contagious disease caused by the variola virus. The disease affected millions of people over thousands of years and variola virus ranked as one of the deadliest viruses in human history. The complete eradication of smallpox in 1980, a major triumph in medicine, was achieved through a global vaccination campaign using a less virulent poxvirus, vaccinia virus. Despite this success, the herd immunity established by this campaign has significantly waned, and concerns are rising about the potential reintroduction of variola virus as a biological weapon or the emergence of zoonotic poxviruses. These fears were further fueled in 2022 by a global outbreak of monkeypox virus (mpox), which spread to over 100 countries, thereby boosting interest in developing new vaccines using molecular approaches. However, poxviruses are complex and creating modern vaccines against them is challenging. This review focuses on the structural biology of the six major neutralization determinants on poxviruses (D8, H3, A27, L1, B5, and A33), the localization of epitopes targeted by neutralizing antibodies, and their application in the development of subunit vaccines.

The poxvirus F17 protein counteracts mitochondrially orchestrated antiviral responses.

Poxviruses are unusual DNA viruses that replicate in the cytoplasm. To do so, they encode approximately 100 immunomodulatory proteins that counteract cytosolic nucleic acid sensors such as cGAMP synthase (cGAS) along with several other antiviral response pathways. Yet most of these immunomodulators are expressed very early in infection while many are variable host range determinants, and significant gaps remain in our understanding of poxvirus sensing and evasion strategies. Here, we show that after infection is established, subsequent progression of the viral lifecycle is sensed through specific changes to mitochondria that coordinate distinct aspects of the antiviral response. Unlike other viruses that cause extensive mitochondrial damage, poxviruses sustain key mitochondrial functions including membrane potential and respiration while reducing reactive oxygen species that drive inflammation. However, poxvirus replication induces mitochondrial hyperfusion that independently controls the release of mitochondrial DNA (mtDNA) to prime nucleic acid sensors and enables an increase in glycolysis that is necessary to support interferon stimulated gene (ISG) production. To counter this, the poxvirus F17 protein localizes to mitochondria and dysregulates mTOR to simultaneously destabilize cGAS and block increases in glycolysis. Our findings reveal how the poxvirus F17 protein disarms specific mitochondrially orchestrated responses to later stages of poxvirus replication.

Genomic characterization and molecular evolution of human monkeypox viruses.

Monkeypox virus is a member of the family Poxviridae, as are variola virus and vaccinia virus. It has a linear double-strand DNA genome approximately 197 kb long, containing ~190 non-overlapping ORFs. Comparison of members of the Central and West African clades shows the presence of unique genes that are associated with different disease presentations, depending on the strain. The last smallpox vaccination efforts ended in the mid-1980s, and there is concern about the recent spread of human monkeypox disease around the world. Almost 87,000 human monkeypox cases have been diagnosed in the world, of which more than 10,900 were in Brazil. The aim of this study was to evaluate the epidemiology and molecular evolution of hMpxV. From computational biology analysis of 640 hMpxV genomes from 1962 to 2022, synteny breaks and gene conservation were observed between Central and West clade genomes, and strains belonged with the 2022 outbreak assigned to the West African clade. Evidence was found for diversifying selective pressure at specific sites within protein coding sequences, acting on immunomodulatory processes. The existence of different sites under diversifying and purifying selection in paralog genes indicates adaptive mechanisms underlying the host-pathogen interaction of monkeypox virus in humans.

Evolution and diversity of nucleotide and dinucleotide composition in poxviruses.

Poxviruses (family Poxviridae) have long dsDNA genomes and infect a wide range of hosts, including insects, birds, reptiles and mammals. These viruses have substantial incidence, prevalence and disease burden in humans and in other animals. Nucleotide and dinucleotide composition, mostly CpG and TpA, have been largely studied in viral genomes because of their evolutionary and functional implications. We analysed here the nucleotide and dinucleotide composition, as well as codon usage bias, of a set of representative poxvirus genomes, with a very diverse host spectrum. After correcting for overall nucleotide composition, entomopoxviruses displayed low overall GC content, no enrichment in TpA and large variation in CpG enrichment, while chordopoxviruses showed large variation in nucleotide composition, no obvious depletion in CpG and a weak trend for TpA depletion in GC-rich genomes. Overall, intergenome variation in dinucleotide composition in poxviruses is largely accounted for by variation in overall genomic GC levels. Nonetheless, using vaccinia virus as a model, we found that genes expressed at the earliest times in infection are more CpG-depleted than genes expressed at later stages. This observation has parallels in betahepesviruses (also large dsDNA viruses) and suggests an antiviral role for the innate immune system (e.g. via the zinc-finger antiviral protein ZAP) in the early phases of poxvirus infection. We also analysed codon usage bias in poxviruses and we observed that it is mostly determined by genomic GC content, and that stratification after host taxonomy does not contribute to explaining codon usage bias diversity. By analysis of within-species diversity, we show that genomic GC content is the result of mutational biases. Poxvirus genomes that encode a DNA ligase are significantly AT-richer than those that do not, suggesting that DNA repair systems shape mutation biases. Our data shed light on the evolution of poxviruses and inform strategies for their genetic manipulation for therapeutic purposes.

Genotyping tool for salmonid gill pox virus (SGPV) obtained from farmed and wild Atlantic salmon (Salmo salar).

Poxviruses are common viruses found in vertebrate species. In 2006, the first poxvirus associated with salmon, salmonid gill poxvirus (SGPV), was identified during an outbreak of gill disease at a smolt production site in northern Norway and at two marine farms in western Norway. Poxviruses had previously been detected in ayu (Plecoglossus altivelis) and koi carp (Cyprinus carpio). In all three fish species, poxviruses are associated with gill disease. It has not been possible to culture SGPV from Norway, and little is known about its virulence. However, the association between SGPV and gill disease in salmon has shown the need for molecular tools to identify reservoirs and transmission routes. Sequencing the genome of a second isolate of SGPV has made it possible to compare variable regions between two strains of the virus, showing the presence of a large number of variable regions that exhibit both variable numbers of tandem repeats and intra-ORF variation. We present eight regions that are suitable for distinguishing strains of SGPV and determining their phylogenetic relationship, and these were used to compare SGPV isolates obtained from both farmed and wild salmon in fresh and sea water. The prevalence of the virus was found to be higher in wild salmon in rivers than in returning wild salmon collected from traps in Norwegian fjords. Genotyping based on the eight selected variable regions, suggests the presence of geographically distinct isolates in freshwater among both farmed and wild salmon, while SGPV from marine farms shows high local diversity and a wide geographical distribution of similar strains of the virus.

Immune responses in carp strains with different susceptibility to carp edema virus disease.

Carp edema virus disease (CEVD), also known as koi sleepy disease (KSD), represents a serious threat to the carp industry. The expression of immune-related genes to CEV infections could lead to the selection of crucial biomarkers of the development of the disease. The expression of a total of eleven immune-related genes encoding cytokines (IL-1ß, IL-10, IL-6a, and TNF-α2), antiviral response (Mx2), cellular receptors (CD4, CD8b1, and GzmA), immunoglobulin (IgM), and genes encoding-mucins was monitored in gills of four differently KSD-susceptible strains of carp (Amur wild carp, Amur Sasan, AS; Ropsha scaly carp, Rop; Prerov scaly carp, PS; and koi) on days 6 and 11 post-infection. Carp strains were infected through two cohabitation infection trials with CEV genogroups I or IIa. The results showed that during the infection with both CEV genogroups, KSD-susceptible koi induced an innate immune response with significant up-regulation (p < 0.05) of IL-1ß, IL-10, IL-6a, and TNF-α2 genes on both 6 and 11 days post-infection (dpi) compared to the fish sampled on day 0. Compared to koi, AS and Rop strains showed up-regulation of IL-6a and TNF-α2 but no other cytokine genes. During the infection with CEV genogroup IIa, Mx2 was significantly up-regulated in all strains and peaked on 6 dpi in AS, PS, and Rop. In koi, it remained high until 11 dpi. With genogroup I infection, Mx2 was up-expressed in koi on 6 dpi and in PS on both 6 and 11 dpi. No significant differences were noticed in selected mucin genes expression measured in gills of any carp strains exposed to both CEV genogroups. During both CEV genogroups infections, the expression levels of most of the genes for T cell response, including CD4, CD8b1, and GzmA were down-regulated in AS and koi at all time points compared to day 0 control. The expression data for the above experimental trials suggest that both CEV genogroups infections in common carp strains lead to activation of the same expression pattern regardless of the fish's susceptibility towards the virus. The expression of the same genes in AS and koi responding to CEV genogroup IIa infection in mucosal tissues such as gill, gut, and skin showed the significant up-regulation of all the cytokine genes in gill and gut tissues from koi carp at 5 dpi. Significant down-regulation of CD4 and GzmA levels were only detected in koi gill on 5 dpi but not in other tissues. AS carp displayed significant up-expression of Mx2 gene in all mucosal tissues on 5 dpi, whereas in koi, it was up-regulated in gill and gut only. In both carp strains, gill harbored a higher virus load on 5 dpi compared to the other tissues. The results showed that resistance to CEV could not be linked with the selected immune responses measured. The up-regulation of mRNA expression of most of the selected immune-related genes in koi gill and gut suggests that CEV induces a more systemic mucosal immune response not restricted to the target tissue of gills.

Structural homology screens reveal host-derived poxvirus protein families impacting inflammasome activity.

Viruses acquire host genes via horizontal transfer and can express them to manipulate host biology during infections. Some homologs retain sequence identity, but evolutionary divergence can obscure host origins. We use structural modeling to compare vaccinia virus proteins with metazoan proteomes. We identify vaccinia A47L as a homolog of gasdermins, the executioners of pyroptosis. An X-ray crystal structure of A47 confirms this homology, and cell-based assays reveal that A47 interferes with caspase function. We also identify vaccinia C1L as the product of a cryptic gene fusion event coupling a Bcl-2-related fold with a pyrin domain. C1 associates with components of the inflammasome, a cytosolic innate immune sensor involved in pyroptosis, yet paradoxically enhances inflammasome activity, suggesting differential modulation during infections. Our findings demonstrate the increasing power of structural homology screens to reveal proteins with unique combinations of domains that viruses capture from host genes and combine in unique ways.

Pandemic potential of poxviruses: From an ancient killer causing smallpox to the surge of monkeypox.

Smallpox caused by the variola virus (VARV) was one of the greatest infectious killers of mankind. Historical records trace back smallpox for at least a millennium while phylogenetic analysis dated the ancestor of VARV circulating in the 20th century into the 19th century. The discrepancy was solved by the detection of distinct VARV sequences first in 17th-century mummies and then in human skeletons dated to the 7th century. The historical records noted marked variability in VARV virulence which scientists tentatively associated with gene losses occurring when broad-host poxviruses narrow their host range to a single host. VARV split from camel and gerbil poxviruses and had no animal reservoir, a prerequisite for its eradication led by WHO. The search for residual pockets of VARV led to the discovery of the monkeypox virus (MPXV); followed by the detection of endemic smallpox-like monkeypox (mpox) disease in Africa. Mpox is caused by less virulent clade 2 MPXV in West Africa and more virulent clade 1 MPXV in Central Africa. Exported clade 2 mpox cases associated with the pet animal trade were observed in 2003 in the USA. In 2022 a world-wide mpox epidemic infecting more than 80,000 people was noted, peaking in August 2022 although waning rapidly. The cases displayed particular epidemiological characteristics affecting nearly exclusively young men having sex with men (MSM). In contrast, mpox in Africa mostly affects children by non-sexual transmission routes possibly from uncharacterized animal reservoirs. While African children show a classical smallpox picture, MSM mpox cases show few mostly anogenital lesions, low-hospitalization rates and 140 fatal cases worldwide. MPXV strains from North America and Europe are closely related, derived from clade 2 African MPXV. Distinct transmission mechanisms are more likely causes for the epidemiological and clinical differences between endemic African cases and the 2022 epidemic cases than viral traits.

Vertebrate-tropism of a cressdnavirus lineage implicated by poxvirus gene capture.

Among cressdnaviruses, only the family Circoviridae is recognized to infect vertebrates, while many others have unknown hosts. Detection of virus-to-host horizontal gene transfer is useful for solving such virus-host relationships. Here, we extend this utility to an unusual case of virus-to-virus horizontal transfer, showing multiple ancient captures of cressdnavirus Rep genes by avipoxviruses-large dsDNA pathogens of birds and other saurians. As gene transfers must have occurred during virus coinfections, saurian hosts were implied for the cressdnavirus donor lineage. Surprisingly, phylogenetic analysis revealed that donors were not members of the vertebrate-infecting Circoviridae, instead belonging to a previously unclassified family that we name Draupnirviridae. While draupnirviruses still circulate today, we show that those in the genus Krikovirus infected saurian vertebrates at least 114 Mya, leaving endogenous viral elements inside snake, lizard, and turtle genomes throughout the Cretaceous Period. Endogenous krikovirus elements in some insect genomes and frequent detection in mosquitoes imply that spillover to vertebrates was arthropod mediated, while ancestral draupnirviruses likely infected protists before their emergence in animals. A modern krikovirus sampled from an avipoxvirus-induced lesion shows that their interaction with poxviruses is ongoing. Captured Rep genes in poxvirus genomes often have inactivated catalytic motifs, yet near-total presence across the Avipoxvirus genus, and evidence of both expression and purifying selection on them suggests currently unknown functions.

ICTV Virus Taxonomy Profile: Poxviridae 2023.

Poxviridae is a family of enveloped, brick-shaped or ovoid viruses. The genome is a linear molecule of dsDNA (128-375 kbp) with covalently closed ends. The family includes the sub-families Entomopoxvirinae, whose members have been found in four orders of insects, and Chordopoxvirinae, whose members are found in mammals, birds, reptiles and fish. Poxviruses are important pathogens in various animals, including humans, and typically result in the formation of lesions, skin nodules, or disseminated rash. Infections can be fatal. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Poxviridae, which is available at ictv.global/report/poxviridae.

Extensive ITR expansion of the 2022 Mpox virus genome through gene duplication and gene loss.

Poxviruses are known to evolve slower than RNA viruses with only 1-2 mutations/genome/year. Rather than single mutations, rearrangements such as gene gain and loss, which have been discussed as a possible driver for host adaption, were described in poxviruses. In 2022 and 2023 the world is being challenged by the largest global outbreak so far of Mpox virus, and the virus seems to have established itself in the human community for an extended period of time. Here, we report five Mpox virus genomes from Germany with extensive gene duplication and loss, leading to the expansion of the ITR regions from 6400 to up to 24,600 bp. We describe duplications of up to 18,200 bp to the opposed genome end, and deletions at the site of insertion of up to 16,900 bp. Deletions and duplications of genes with functions of supposed immune modulation, virulence and host adaption as B19R, B21R, B22R and D10L are described. In summary, we highlight the need for monitoring rearrangements of the Mpox virus genome rather than for monitoring single mutations only.

Detection of Novel Poxvirus from Gray Seal (Halichoerus grypus), Germany.

We detected a novel poxvirus from a gray seal (Halichoerus grypus) from the North Sea, Germany. The juvenile animal showed pox-like lesions and deteriorating overall health condition and was finally euthanized. Histology, electron microscopy, sequencing, and PCR confirmed a previously undescribed poxvirus of the Chordopoxvirinae subfamily, tentatively named Wadden Sea poxvirus.

Improving poxvirus-mediated antitumor immune responses by deleting viral cGAMP-specific nuclease.

cGAMP-specific nucleases (poxins) are a recently described family of proteins dedicated to obstructing cyclic GMP-AMP synthase signaling (cGAS), an important sensor triggered by cytoplasmic viral replication that activates type I interferon (IFN) production. The B2R gene of vaccinia viruses (VACV) codes for one of these nucleases. Here, we evaluated the effects of inactivating the VACV B2 nuclease in the context of an oncolytic VACV. VACV are widely used as anti-cancer vectors due to their capacity to activate immune responses directed against tumor antigens. We aimed to elicit robust antitumor immunity by preventing viral inactivation of the cGAS/STING/IRF3 pathway after infection of cancer cells. Activation of such a pathway is associated with a dominant T helper 1 (Th1) cell differentiation of the response, which benefits antitumor outcomes. Deletion of the B2R gene resulted in enhanced IRF3 phosphorylation and type I IFN expression after infection of tumor cells, while effective VACV replication remained unimpaired, both in vitro and in vivo. In syngeneic mouse tumor models, the absence of the VACV cGAMP-specific nuclease translated into improved antitumor activity, which was associated with antitumor immunity directed against tumor epitopes.

Exaptation of Inactivated Host Enzymes for Structural Roles in Orthopoxviruses and Novel Folds of Virus Proteins Revealed by Protein Structure Modeling.

Viruses with large, double-stranded DNA genomes captured the majority of their genes from their hosts at different stages of evolution. The origins of many virus genes are readily detected through significant sequence similarity with cellular homologs. In particular, this is the case for virus enzymes, such as DNA and RNA polymerases or nucleotide kinases, that retain their catalytic activity after capture by an ancestral virus. However, a large fraction of virus genes have no readily detectable cellular homologs, meaning that their origins remain enigmatic. We explored the potential origins of such proteins that are encoded in the genomes of orthopoxviruses, a thoroughly studied virus genus that includes major human pathogens. To this end, we used AlphaFold2 to predict the structures of all 214 proteins that are encoded by orthopoxviruses. Among the proteins of unknown provenance, structure prediction yielded clear indications of origin for 14 of them and validated several inferences that were previously made via sequence analysis. A notable emerging trend is the exaptation of enzymes from cellular organisms for nonenzymatic, structural roles in virus reproduction that is accompanied by the disruption of catalytic sites and by an overall drastic divergence that precludes homology detection at the sequence level. Among the 16 orthopoxvirus proteins that were found to be inactivated enzyme derivatives are the poxvirus replication processivity factor A20, which is an inactivated NAD-dependent DNA ligase; the major core protein A3, which is an inactivated deubiquitinase; F11, which is an inactivated prolyl hydroxylase; and more similar cases. For nearly one-third of the orthopoxvirus virion proteins, no significantly similar structures were identified, suggesting exaptation with subsequent major structural rearrangement that yielded unique protein folds. IMPORTANCE Protein structures are more strongly conserved in evolution than are amino acid sequences. Comparative structural analysis is particularly important for inferring the origins of viral proteins that typically evolve at high rates. We used a powerful protein structure modeling method, namely, AlphaFold2, to model the structures of all orthopoxvirus proteins and compared them to all available protein structures. Multiple cases of recruitment of host enzymes for structural roles in viruses, accompanied by the disruption of catalytic sites, were discovered. However, many viral proteins appear to have evolved unique structural folds.

Characteristics of pathology and transcriptome profiling reveal features of immune response of acutely infected and asymptomatic infected of carp edema virus in Koi.

Koi sleepy disease (KSD) is a high mortality and infection viral disease caused by carp edema virus (CEV), which was a serious threat to aquaculture of common carp and export trade of Koi worldwide. Asymptomatic infection is an important cause of the difficulty in preventing KSD and its worldwide spread, because asymptomatic infection can be activated under appropriate condition. However, the understanding of the molecular correlates of these infections is still unknown. The purpose of this study was to compare the pathology change, enzyme activity, immunoglobulin activity, host and viral gene expression differences in acutely infected and cohabiting asymptomatic Koi infected with CEV. Healthy Koi were used as a control. The gross pathology, histopathology and ultrastructural pathology showed the difference and characteristics damage to the tissues of Koi under different infection conditions. Periodic Acid-Schiff stain (PAS), enzyme activity and immunoglobulin activity revealed changes in the immune response of gill tissue between acutely infected, asymptomatic infected and healthy Koi. A total of 111 and 2484 upregulated genes and 257 and 4940 downregulated genes were founded in healthy Koi vs asymptomatic infected Koi and healthy Koi vs acutely infected Koi, respectively. Additionally, 878 upregulated genes and 1089 downregulated genes were identified in asymptomatic vs. acutely infected Koi. Immune gene categories and their corresponding genes in different comparison groups were revealed. A total of 3, 59 and 28 immune-related genes were identified in the group of healthy Koi vs asymptomatic infected Koi, healthy Koi vs acutely infected Koi and asymptomatic infected Koi vs acutely infected Koi, respectively. Nineteen immune-related genes have the same expression manner both in healthy Koi vs acutely infected Koi and asymptomatic Koi vs acutely infected Koi, while 9 immune-related genes were differentially expressed only in asymptomatic Koi vs acutely infected Koi, which may play a role in viral reactivation. In addition, 8 differentially expressed genes (DEGs) were validated by quantitative reverse transcription PCR (RT-qPCR), and the results were consistent with the RNA-Seq results. In conclusion, the data obtained in this study provide new evidence for further elucidating CEV-host interactions and the CEV infection mechanism and will facilitate the implementation of integrated strategies for controlling CEV infection and spread.